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ABSTRACT 



The dark matter (DM) halos of field elliptical galaxies have not been well studied and their properties appear controversial in the 
literature. While some galaxies appear to be nearly devoid of DM, others show clear evidence of its presence. Furthermore, modified 
Netonian dynamics (MOND), which has been found to have predictive power in the domain of disk galaxies, has not yet been 
investigated for isolated elliptical galaxies. We study the kinematics of the isolated elliptical NGC 7507, which has been claimed 
as a clear example of DM presence in early-type galaxies. We obtained major and minor axis long-slit spectroscopy of NGC 7507 
using the Gemini South telescope and deep imaging in Kron-Cousins R and Washington C using the CTIO/MOSAIC camera. Mean 
velocities, velocity dispersion and higher order moments of the velocity distribution are measured out to ~90". The galaxy, although 
almost circular, has significant rotation along the minor axis and a rapidly declining velocity dispersion along both axes. The velocity 
dispersion profile is modeled in the context of a spherical Jeans analysis. Models without DM provide an excellent representation 
of the data with a mass-to-light ratio (M/L) of 3.1 (R-band). The most massive Navarro-Frenk- White (NFW) halo the data allow 
has a virial mass of only 3.9^'j x 10" Mq, although the data are more consistent with models that have a slight radial anisotropy, 
which implies the galaxy has an even lower DM halo mass of 2.2+ j'? x 10" Mq. Modeling of the /14 Gauss-Hermite coefficient 
is inconclusive but seems to be consistent with mild radial anisotropy. A cored logarithmic DM halo with parameters r = 7 kpc 
and Vo = 100 km s can also reproduce the observed velocity dispersion profile. The MOND predictions overestimate the velocity 
dispersion. In conclusion, we cannot easily reproduce the previous findings of a predominance of DM in NGC 7507 within a simple 
spherical model. DM may be present, but only in conjunction with a strong radial anisotropy, for which there are some indications. 

Key words. Galaxies: individual: NGC 7507 - Galaxies: kinematics and dynamics 



1. Introduction 

The existence of DM in spiral galaxies (or alternatively a law of 
gravitation that differs from the Newtonian law) remains undis- 
puted. However, the same cannot be said conclusively for all 
elliptical g alaxies. Central ga laxies in gal axy clusters such as 
NGC 6166 dKelson et al. 2002]), NGC 1 399 dRichtler et alj2008t 
ISchuberth et al.ll2010l) . and NGC 3311 dRichtler etalJl20 lib are 
embedded in massive dark halos as is the case of bright non- 
central cluster ellipticals such as NGC 4636 and NGC 4374 in 
Fornax (Schu berth et al.ll2006t iNapolitano et al.ll20ldl) . In con- 
strast elliptical galaxies with lower luminosities in sparser en- 
vironments (NGC 821, 3379, 4494, 4697), whose mass profiles 
have been studied using planetary nebulae ( PNe) kinematics, ap- 
pear to show only weak evidence of DM (Romanowsky et al] 



2003t [Douglas et al.ll2007HNapolitano etaill2009HMendez et all 



* Based on observations taken at the Gemini Observatory, which is 
operated by the Association of Universities for Research in Astronomy, 
Inc., under a cooperative agreement with the NSF on behalf of the 
Gemini partnership: the National Science Foundation (United States), 
the Science and Technology Facilities Council (United Kingdom), 
the National Research Council (Canada), CONICYT (Chile), the 
Australian Research Council (Australia), Ministerio da Ciencia e 
Tecnologia (Brazil) and SECYT (Argentina). 

** Current address: Finnish Centre for Astronomy with ESO 
(FINCA), University of Turku, Vaisalantie 20, FI-21500 Piikkio, 
Finland; rcsave@utu.fi 



2009; iTeodorescu et al.ll2010l but see lWeiimans et al.ll2009l) . al- 
though the case of NGC 3379 may be consistent with a nor- 
mal cold DM halo if a very strong radial anisotropy is present 
dde Lorenzi et alj|2009l) . 

The comparison of field, and possibly isolated, elliptical 
galaxies to those in clusters is of great interest. Different for- 
mation mechanisms may have led to differences in their dark 
halo structure, but observationally few properties have been re- 
liably determined. Isolated galaxies may also be goo d test ob- 
jects of the Modified Newtonian Dynamics (MOND, Milgrom 
11983b ISanders & McGaugrill2002l). since complications arising 
from the "external field effect" dMilgromlll983t IWu et alj|20071) 
should be minimal. 

The galaxy NGC 7507 is an elliptical (E0, Ide Souza et ail 
120041) of M y =-21.6, with a distance modulus of 31.83+0.17 
based on s urface brightness fluctuations (23.22±1.8 Mpc; 
i Tonrv et ail 1200 ll corrected by -0.16 mag following 
iJensen et al.ll2003l) ; we assume this distance throughout the pa- 
per, which implies a scale of 112.5 pc arcsec -1 . The galaxy is 
fairly isolated, forming a probable pair with the SBb galaxy 
NGC 7513, which lies at a projected distance of 18 arcmin- 
utes ( ~ 130 kpc), but displaying no signs of recent interac- 
tions dTal et alJl2009l) . Due to its proximity and regular shape, 
its dynamics have been well studied l eading to strikingly dif- 
ferent conclusions. Berl in et alJ d 19941) found no evidence of 
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a DM halo, using maj or-axis kinematical data out to 66". 
iRronawitter et alJ (|2000, hereafter K+00) performed dynami- 
cal analyses of 21 round elliptical galaxies, including NGC 
7507. They concluded that models based on luminous mat- 
ter could only be ruled out with 95% confidence for only 
thre e galaxies, including NGC 75 07. Contradicting this find- 
ing, Magor rian & Ballantvnd (1200 ll) adopted the kinema tics of 
Berti n et alJd!994t) and photometry of lFranx etall d!989bl) to es- 
tablish that NGC 7507 has a constant mass-to-light ratio (MIL) 
profile, i.e. that it does not contain a dark halo. 

The most significant difference between these studies ap- 
pears to be the observational data used and not the modeling 
approach, which, although different, should not give strongly 
discrepant results. We decided to obtain new, deeper long-slit 
spectroscopy using Gemini South. Eight-meter class telescopes 
have been scarcely used to measure the long-slit kinematics 
of the unresolved stellar components of elliptical galaxies (e.g . 
iThomas et al]l2007t iForestell & Gebhardill2()l(t iPu et aIll2010L 
and see Samurovic & Danzigerll2005l for a review of the work 
made with 4m class telescope s). Much effort ha s been put into 
kinematic studi es of PNe (e.g.lMendez et al. 2009) an d globula r 
clusters (e.g. iRomanowskv et all 120091: ISchuberth et al J 12010), 
which can be observed out to much larger galactocentric radii, 
but provide weaker constraints of their orbital anisotropies. 

The paper is organized as follows: we present our observa- 
tions and reduction procedures in Sect. [2] The photometric prop- 
erties of the galaxy together with a description of its kinemat- 
ics are given in Sect. [3] The dynamical modeling is presented 
in Sect. |U while the discussion and conclusions can be seen in 
Sects. |5]and|6j respectively. 



2. Observations and data reduction 

2.1. CTIO/MOSAIC photometry and photometric calibration 

The images used in this study were obtained with the MOSAIC 
II camera (an 8 CCDs mosaic imager, 16 amplifiers) mounted 
at the prime focus of the 4-m Blanco telescope at the Cerro 
Tololo Inter- American Observatory (CTIO, Chile). Observations 
were carried out during 2005 August 5-6, where the first night 
was "useful" and the second one "photometric", according to 
the CTIO reports. One pixel of the MOSAIC wide-field camera 
subtends 0.27 arcsec on the sky, which corresponds to a field of 
36'x 36', i.e. approximately 230 x 230 kpc 2 at the distance of 
NGC 7507. 

All fields were imaged in Washington C (/t cen trai = 3850 A) 
and Kron-Cousins R (/i ce ntrai = 6440 A) filters. The R fil t er was 
used instead of the original Washington T\, as Geisler (1996) 
has shown that the Kron-Cousins R filter is more efficient than 
T\, and that R and 7T magnitudes are closely related, with just a 
very small colour term and zero-point difference (R—Tl — 0.02). 

To fill in the gaps between the eight individual MOSAIC 
chips, the data were dithered taking four images in R with expo- 
sure times of 720 s each, and seven images in C with exposures 
of 1800 s each. To avoid saturation at the galaxy core, additional 
shorter exposures of 60 s in R and of 300 s in C were also ob- 
tained. Sky flats were obtained at the beginning and end of each 
night. 

The MOSAIC data reduction was performed using tasks 
from the MSCRED package within IRAF. Several tasks were em- 
ployed to correct for the variable pixel scale across the CCD, 
which might cause a 4 per cent variability in the point-sources 
brightness (from the center to the corners). The final combined 
R image still shows sensitivity variations up to 1 .4 per cent, and 



the final C image up to 3.0 per cent (peak-to-peak). The seeing 
in these final images is 1.1" on the R frame, and 1.3 " on the C 
image. 

The photometric calibration was perfo rmed with three fields 
of standard stars selected from the list of Geisler (1996), each 
containing 7 to 10 standard stars. For these fields, C and R im- 
ages were obtained in each night, covering an airmass range 1 . 1 
- 2.3. Finally, we used a single set of calibration equations for 
both nights. 

The equations read 

Tl = R inst + (0.89 ± 0.01) - (0.23 ± 0M)X R 
+(0.016 + 0.003)(C - 7T), 
C = C inst + (1.55 ± 0.01) - (0.30 + 0.01)X C 
+(0.069 + 0.003)(C - 7T), 

where C; nst and /?i nst are instrumental magnitudes and Xq and Xr 
are the airmasses. 

The absolute photometric quality of the C-photometry is du- 
bious, hence it is only used for differential analyses. The R- 
photometry is of high quality, as the comparison of our results 
with previous photomet ric data se t s of the galaxy reveals. The 
mean difference from iFranx et a D d~L989b) is 0.08 mag (Franx 
et al. is brighter). We a lso compare with photo electric aperture 
measurements given bv lPoulain & Nietol (l994). In this case, we 
find that our profile is brighter by 0. 15 mag. 

2.2. Surface photometry 

The surface photometry of the galaxy was perfomed using the 
IRAF/ELLIPSE task. The shallower images were used to deter- 
mine the surface brightness profile for the inner ~10" in the C 
and R profiles, which were saturated in the longer exposures. The 
sky brightness in the shallow images was determined by taking 
the median value of the sky in several 15 x 15 pixel boxes of 
empty regions further than 5' from the galaxy center. Since an 
accurate sky determination is critical in shaping the light pro- 
file in the outer parts of a galaxy, a differ ent procedure, similar 
to the one devised by iDirsch et al.1 d2003 | ). w as adopted for the 
deep images. First, DAOPHOT 7stetsonlll98"7h was applied to all 
point sources present in the field. The local sky brightness of 
each source was determined using the median of the sky value 
of the pixels between 15 and 25 pixels from each source center. 
Deviant pixels were rejected using a three-sigma clipping pro- 
cedure. All the sky values were then ordered according to their 
distance from the center of NGC 7507. To exclude overestimated 
sky values (of point sources close to nearby galaxies, bright stars 
or chip defects), a new iterative 3-sigma clipping was applied to 
the sky values within radial bins of 100 objects, t his time us- 
ing ro bust location and scale bi-weight estimators dBeers et al.l 
1990). At distances larger than 5' from the center of the galaxy, 
the value of the sky brightness fluctuates by less than 0.5%. The 
final value for the sky was then taken as the mean of the sky val- 
ues between 10' and 12' (~78-98 kpc) from the galaxy center 
and the R sky brightness in mag/arcsec 2 was found to be 20.90. 
The photometry of the galaxy can be seen in Table Q] 

2.3. Gemini/GMOS spectroscopy 

Long-slit spectra of NGC 7507 were obtained using the Gemini 
Multi-Object Spectrograph (GMOS) on Gemini South, Cerro 
Pachdn, Chile, in queue mode during the nights of September 
19 and 21, 2004 (program GS-2004B-Q-75). The grism was the 
B600+_G5323 which, in conjunction with a 1" slit width, gives 
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a resolution of ~4.7A FWHM. Three consecutive exposures of 
1800s each were taken for a position angle of 90 °, which is 
close to the major axis (PA = 105°), and for PA - 0°, which 
is close to the minor axis (PA = 15°). For convenience, these 
two positions are referred to hereafter as the major and minor 
axes, respectively. Standard reduction steps (bias subtraction, flat 
fielding, detectors mosaicing, wavelength calibration) were per- 
formed with the Gemini IRAF package (v 1.9.1). Typical wave- 
length residuals were ~0.05 A. No cosmic ray cleaning was at- 
tempted at this stage because it introduced undesired noise, es- 
pecially around bright sky lines. Since copper-argon arc spectra 
were taken during daytime and not immediately either after or 
before the science exposures, a correction of the zero point of 
the wavelength solution was necessary and was performed by 
measuring the position of the bright sky line at 5577.34 A. This 
correction typically amounted to 0. 1 A. After that, the three spec- 
tra for each position angle were averaged. Sky subtraction was 
done using the IRAF task background, by averaging 50 rows 
on opposite sides of each spectra, at around 4.5' from the galaxy 
center for the major axis and ~ 4' for the minor axis. The two- 
dimensional (2D) spectra were spatially rebinned to achieve a 
constant S /N ~ 50/pixel, which is highly desirable when mea- 
suring the higher order moments of the line-of-sight velocity dis- 
tribution (LOSVD). The 4800-5500 A wavelength range, which 
contains strong absorption features as Hfi, Mg b, and Fe lines 
around 5325 A, was selected for the kinematical analysis. 

Mean velocity, velocity dispersion and higher order moment 
profiles were measured for both long slits by u sing the Gauss- 
Hermite pixel fitting code (Ivan der MarellH994l hereafter vdM 
co de) and the penalized pixel fittin g (pPXF) method developed 
by ICappellari & Emsellem (2004). Both approaches have the 
advantage of performing their calculations in pixel space in- 
stead of Fourier space, thus allowing a more precise masking 
of undesired zones in the spectra (emission lines, bad columns, 
etc). Both methods also parametrize the LOSVD in terms of 
a Gauss-Hermite series dGerhardlll99"l l van der Marel & Franxl 
1 1993b . The former algorithm measures the mean velocity and 
the velocity dispersion, and, with these values fixed, fits the 
higher order terms of the velocity distribution, parametrized by 
the Gauss-Hermite coefficients hj, and h\. The latter algorithm 
fits all parameters simultaneously, but can bias the solution to a 
Gaussian shape when th e higher order terms are not s ufficiently 
constrained by the data (Cappellari & Emsellem 2004). An ad- 
vantage of the pPXF method over the vdM code is that it can 
perform the fitting process with a linear combination of template 
spectra, minimizing the effect of template mismatching. 

A subset of 22 old and metal-r ich single stellar population 
models from Vazdekis et al. (2010) based on the MILES stel- 
lar library (Sa nchez-Blazquez et al.ll2006l) were used as velocity 
templates for our analysis using the vdM code. These templates 
were logarithmically rebinned to the same dispersion and con- 
volved with a Gaussian to ensure that they had the same spectral 
resolution as the science spectra. Final values of mean velocity, 
velocity dispersion, hs and h\ parameters were assumed to be 
the average of the values for the three templates that provided 
the closest fit (the lowest^ 2 values). 

In the pPXF analysis, two different sets of spectra were em- 
ployed as templates to extract the kinematical information, the 
aforementioned Vazdekis et al. models, but this time the en- 
tire library c onsisting of 350 models , and t he full MILES spec- 
tral library dSanchez-Blazquez et alJ 12006). Both have a spec- 
tral resolution of 2.3 A FWHM (but see iBeifiori et alJ 1201 ll) . 
Since we expect to find a population gradient across the galaxy 
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Fig. 1. Comparison between the different measurements of the 
velocity dispersion profile of the major axis of NGC 7507 as 
described in the text. 



(see Section O, the templates used for each radial bin were not 
fixed, but freely chosen by pPXF from the MILES library and 
the Vazdekis models, separately, until the closest model was 
achieved. For the Vazdekis models, pPXF used typically 3 from 
the 350 models, while a mean number of 16 template stars from 
the 985 in the MILES library were used by pPXF for each ra- 
dial bin. Estimations of the errors for all the calculated values 
were obtained by performing 100 Monte Carlo simulations in 
which random noise was added to eac h science spectrum and 
then passed again th rough pPXF (e.g. Cappellari & Emsellem 
12004 iBedregal et aljr2006l) . In this case, only the templates that 
were previously selected by pPXF were fitted, and not the entire 
set. The dispersion in the values obtained from the simulations 
are quoted as the error bars in the measurements from the origi- 
nal spectra. 

The agreement between the results of pPXF using the two 
sets of templates is very good. However a few remarks are ap- 
propriate. Typically pPXF uses much fewer templates from the 
Vazdekis models (~3) and the x 1 values are always larger than 
when using MILES. The inclusion of multiplicative Legendre 
polynomials reduces the y 2 , but does not increase the number 
of templates used. We find a systematic small difference of 3 
km s" 1 between the velocity dispersion results obtained using 
the Vazdekis SSP models and the MILES library (Fig. Q] bot- 
tom panel). The difference from the values given by the vdM 
code is somewhat larger: the velocity dispersions derived with 
vdM are systematically larger. We find a systematic difference 
of 9 km s _1 (Fig. Q] top panel). Since template mismatching is 
expected to be far less significant for this method, if not totally 
negligible, based on p revious studies (e.g. IShapiro et alj [2006: 
Cappellari et al. 2007), we decided to adopt the results obtained 
using the MILES spectral in all the subsequent kinematic de- 
scription and dynamical analyses. These measurements can be 
seen in Tables [2] and [3] 
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Fig. 2. CTIO/MOSAIC surface photometry of NGC 7507. a) 
Ellipticity profile, b) C — R color profile, c) R surface bright- 
ness profile (black circles) and R surface photometry from 
iFranx et all d 198 9b) (filled red squares). The solid red line is the 
fit of a double B model to our data. Residuals can be seen in 
panel d. 



3. Results 

3. 1 . Light and color profile 

There is a well-established variety of analytical model light pro- 
files of elliptical galaxies. Since we wish to use the light pro- 
file within a spherical Jeans analysis, its analytical form, besides 
being an excellent representation, should permit a closed ana- 
lytical deprojection as well as an analytical mass profile. These 
conditions are fullfilled by the sum of two B profiles (or Hubble- 
Reynolds profile) 



H(R) = -2.5 log^a 



l + l* 

Ri 



2\ 



+ a 2 



-0, 



(i) 



where a, = 2.052xl(T 8 , a 2 = 6.452xl(T 7 , Ri = 16", and 
/?2=2.19", provide an excellent representation of the surface 
brightness profile between 2"and 400"(Fig.[2] panel c). 

Although the j6-profile provides an excellent description out 
to about 50 kpc, the extrapolation to larger radii, which is used in 
the dynamic modeling, is not a priori justified. It is generally un- 
known where a galaxy "ends" (but see the remarks in Sect. 14. II ). 
The mass of the /J-model diverges for large radii. This raises the 
concern of whether systematic errors are introduced by extrapo- 
lating the profile to larger radii. 

We therefore fit the surface brightness profile with the sum 
of two Sersic profiles 

I(R) = I exp(-(R/a s ) 1/m ) (2) 

(ISersidfT9 68: Graham & Driverll2 005). since a single Sersic pro- 
file fails to accurately reproduce the luminosity in the entire 
radial range, underestimating the outer (> 20 kpc) profile. The 




Fig. 3. Surface brightness profile of NGC 7507. The green dots 
indicate the R photometry in units of L Q pc~ 2 . The dashed blue 
lines represent the two Sersic profiles, while the solid blue line 
is their sum. The red solid line represents the double B model. 



Sersic parameters are I' Q = 1.90 x 10 6 L pc 2 , a' s = 0.0677 pc, 
to' = 4.8, for the inner Sersic profile, and IS. = 15.77 L pc -2 , 
a° = 15832 pc, m° = 1.05 , for the outer profile; M 0R =4.42 
dBinnev & M errifield fl998l) was adopted. 

Both models provide excellent descriptions of the surface 
density profile from ~ 200 pc (~1.5") until the last photometric 
point at ~ 55 kpc, and are practically indistinguishable (see Fig. 

Differences are seen only in the inner ~200 pc where the dou- 
ble B model underestimates the light profile, while the opposite 
happens with the double Sersic model and for the extrapolation 
of the light profile (outside ~ 55 kpc) where the double B model 
falls as R~ 2 , but the Sersic profile declines more steeply. 

The double Sersic profile implies that the effective radius is 
R e -15 " . This is significantly larger than t he previous estima tions 
of 24"(Lauberts & Valentijn 1989), 31" dFaber et al.ll 19891) . and 
47" dde Souza et aLll2004l) ; this is unsurprising given our deep 
imagi ng. For example, the value of Faber et al. is taken from the 
RC2 dde Vaucouleurs et al.l[l9 76) where the total galaxy bright- 
ness is defined as the brightness within an isophote with a sur- 
face brightness in the fi-band of 25 mag arcsec -2 which we 
found for NGC 7507 occurs at a distance of about 150". 

The color profile was evaluated as the difference between 
the two models in C and R (Fig. [2] panel b). The inner core 
(within 10") declines by 0.3 mag with radius to bluer colors, 
followed by a radial interval of more or less constant color. 
The color then gets redder again at about 50", and remains 
practically unchanged out to 200". The outer gradient to bluer 
colors might exist, but owing to the large uncertainties cannot 
be claim ed convincingly. Co lor p rofiles for the inner 50" were 
given bv lFranxetall (Il989bl) and iGoudfrooii et alj d!994l) . The 
inner strong gradient is apparent in both of them, but the gradi- 
ent of almost 0.5 mag in B - I between 30" and 60" seen by 
IGoudfrooii et al.1 (1 19941) is not found in our color profile. 

It is interesting to consider a 2D color image of the inner 
region. Fig. |4]was obtained by dividing the short C by the R ex- 
posures. It is apparent that the "color field" is not spherically 
symmetric but tends to have bluer colors along the E-W axis. 
One could therefore interpret the plateau with C-R-2.Q between 
10" and 50" as a result of averaging along isophotes, before the 
color gets redder again. We are unable to determine, however, 
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Fig. 4. C - R color image of NGC 7507 based on the short ex- 
posures. Image has been median smoothed in 3 x 3 pixel boxes. 
Image size is 2' x 2'. 



whether diffuse dust or population effects or both are responsi- 
ble. 



3.2. Ellipticity and higher isophote moments 

As all previous studies have pointed out, the galaxy in projection 
is very round. We find (e) = 0.045 + 0.008, where the mean is 
taken over the inner 2', avoiding the innermost 2" which are 
affected by seeing. Beyond 2' the ellipticity slightly increases, 
partly due to the difficulty in masking a bright 11 th mag star, 
located ~6' SE of the galaxy (Fig|2] panel a). 

The Fourier coefficient 04 describes the isophote's boxiness 
or diskyness (negative: boxy; positive:disky). It is zero out to 2'. 
The slight boxiness at larger radii may also be caused by the 
bright star. A table with the photometry, including ellipticities, 
PA, and 04 is available online. 



3.3. Irregular features 

A small, slightly off-centered dip in the light can be seen in 
the shallow R image. It is not apparent in the shallow C image, 
which has poorer seeing. We i nterpret this as a dust lane already 
det ected bylFranx et al.Tdl989bl) . 

iTal et al] (|2009) found evidence of a faint shell north of the 
galaxy. However, when the smoothed 2D model created with 
IRAF\BMODEL task is subtracted from the galaxy light, we see 
thi shell in neither the C nor in the R images. No substruc- 
ture was detected e ither in the shallower images analyzed by 
IdeSouza etaTI(l2004l) . 



3.4. Stellar kinematics 

The kinematics along the majo r axis of NGC 7507 have been 
measured by several groups dBertin et al.l Il994t iFranx et aD 



Fig. 5. Major axis kinematics of NGC 7507, as extracted with 
pPXF (see text), a) Hermite coefficient h\ profile, b) /13 profile, 
both panels show a dashed line at the zero level for compari- 
son, c) velocity dispersion. The horizontal error bars show the 
regions across which the spatial binning was done (shown only 
in this panel, but valid for all four). The vertical dashed strip in- 
dicates a gap between the detectors in which no measurements 
are possible, d) rotation profile. For all panels, open black circles 
are our d ata, red solid square s are from lK+OCX. solid green circles 
are from Franx et al. ( 1989a), and blue open squares come from 
iBertin et al.l(11994l) . 



?9al K+00), but along the minor axis only b y IFranx et aD 
'891 ; 



1989a 

( 1989a) and higher or der moment s only b v lK+OOt We extended 
the radial coverage of IFranx et al .1 (Il989al) by a factor of 1 .5 for 
the major axis and 2.5 for the minor axis. 



3.4.1. Rotation 

The major axis rotation is very low (Fig. [5] panel d), consis- 
tent with no rotation, with a formal amplitude of -1 ±6 km s _1 
within the inner 30". In the 30-60" range, the amplitude slightly 
increases to a value of -13 + 7 km s _1 , which is consistent with 
the aforementioned previous studies. The outermost measured 
bin, at a mean distance of 97"shows a significant increase in the 
rotation with respect to the center (-68 + 8 km s _1 ). This may be 
partly due to the misalignment of the slit direction and the pho- 
tometric major axis but perhaps also a hint that the kinematic 
properties of the galaxy suffer a change at large radii, as se en for 
example in NGC 821 and NGC 1400 dProctor et al.ll2009h . 

The minor axis rotation is pronounced and slightly asym- 
metric with amplitudes of 30 ± 7 km s _1 in the west direction 
and -43 + 7 km s -1 in the east directio n inside the range 5- 
50"(Fig. |6] panel d). IFranx et all {l989a) measured an ampli- 
tude of 36 ± 5 km s _1 in the range 5-30", without detecting any 
asymmetry. In the outermost east and west bins, around 60", 
the mean velocity is consistent with no rotation. This rotation 
profile closely resembles the major axis rotation of NGC 3379 
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Fig. 6. Minor axis kinematics of NGC 7507. The vertical dashed 
strips show the positions of the gaps between detectors and two 
interloping background galaxies, which were masked in the spa- 
tial binning process. Symbols are the same as in Fig. [5] 



dStatler & Smecker-Handfl999l; IWeiimans etai1l2009l their Fig. 
4). 

Summarising, the isophote shapes are sufficiently undis- 
turbed and the projected rotation amplitude is small enough to 
justify the use of a non-rotating, spherical mass model. 

3.4.2. Velocity dispersion 

The major axis velocity dispersion rapidly declines in the inner 
30"and undergoes a slower decline at larger radii. The disper- 
sion values of K+O Oj deviate stron gly, not only from our v alues, 
but also from lFranx etail (Il989al) and iBertin et"ai1 d 1994 (Fig. 
[5] panel c). While the difference at the center of the galaxy is 
roughly within errors (9.7 + 4.8 km s~'), and it can be easily 
explained by the spectra being taken under different seeing con- 
ditions; at 27" from the c enter, th e difference amounts to 70 ± 
14 km s _1 in the sense that iK+OOl values are larger than ours. 

The minor axis velocity dispersion profile shows a similar 
rapid decline in the inner 30" but a less steep decline at larger 
radii (Fig. [6] p anel c). As in th e case of the major axis, the differ- 
ence from lBertin et al.l(ll994l) is minimal and amounts to 14 + 24 
km s" 1 , where the uncertainty is caused almost entirely by the 
IBertin et all d!994h values. 

A kinematic classification of early-type gala xies, based on 
an ang ular momentum proxy was introduced bv lEmsellem et all 
(2007). They defined 



2 FiRi \Vi\ 



(3) 



V 2 



where F, is flux, Rt is distance to the center of the galaxy, 
and Vj and <x,; are the mean velocity and velocity dispersion. 
Emselle m et al.l d2007l) found that early-type galaxies can be sep- 
arated into two classes of slow and fast rotators, depending on 



whether they have Ar values below (slow) or above (fast) 0.1 
inside R e . Even though this parameter was introduced for the 
S AURON 2D data, it can been applied to long-slit measurements 
in order to help us to compare with the o riginal definition, once 
a proper scaling fac tor of 0.57 is applied ( Cappell ari et al . 2007; 
ICoccato et"aT] |2009). In that case the sum is over the quantities 
along the slit. 

We calculated the Ar parameter separately for each long slit 
position. For both axes, the parameter mildly increases inside 
R e reaching A™ n = 0.088 and A™> = 0.018, staying basically flat 
beyond that radius. This is consistent with the mild rotation 
of the galaxy. NGC 7507 can be classified as a slow rotator, 
sharing also the kinematic misalignment which is often seen 
in these galaxies, although it has a rapidly declining velocity 
dispersion profile wh ich is more common amongst fast rotators 
dCoccato et al.ll2009l) . 

3.4.3. Higher order moments 

The Gauss-Hermite coefficient which measures the asym- 
metric deviation from a Gaussian, is consistent with being zero 
along the major axis. Along the minor axis within 20", it appears 
to have a gradient which would be consistent with its rotation 
(Fig. [6] panels b and d). 

The coefficient h\, which indicates symmetric deviations 
from a Gaussian profile, is positive for both the major an d minor 
axis spectra. The sign agrees with the measurements of IK+OOl 
who also found slightly positive values of h\, although their val- 
ues were somewhat smaller (Fig. [5] panel a). We return to this 
result when discussing any possible anisotropies. 



4. Dynamical modeling 

4.1. The stellar mass profile 

To obtain the (spherical) stellar mass profile we have to depro- 
ject the surface brightn ess profile after corr ecting for Galactic 
absorption (A R = 0.128; ISchlegel et al.|[l998l) . In the case of the 
double-/? profile this is achieved by applying the Abel deprojec- 
tion formula, which for the ft profile has an analytical description 



dt/Ri 



[l + ir/Rd 2 ] 3 ' 2 ' 



(4) 



where a, and R, come from Eq.Q] and Cr = 2.5 x 10 10 is the fac- 
tor used to transform to uni ts of L.0 pc~ 2 , adopting Mq r = 4.42 
dBinnev & Merrifieklll998h . The enclosed stellar mass can be 
obtained by integrating this light density 



M*(r) = An f Y*,«(r);(r)r 2 dr, 

Jo 



(5) 



where Y+.« is the stellar M/Lr. For constant M/Lr, this integral 
has the exact solution 



M. 



.(r) = 27rC R Y^^a,/? 2 



sinh Xj - 



Xi 



ft 



+ 1 



(6) 



where Xt = r/Rj and a,- and R, have the same meaning as in Eq. 

El 

The stellar M/Lr appears as a parameter in the fit to the kine- 
matical data, adopting the deprojected light profile and certain 
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forms for a possible DM halo. From the shape of the colour pro- 
file, one concludes that M/L as a population property is approx- 
imately constant out to a radius of 2'. At larger radii, the M/L 
may drop owing to a declining metallicity. However, this is only 
weakly constrained and any effect on the inner projected velocity 
dispersion would be completely degenerate with a dark halo. 

The Sersic profile can be deproi ected using an approximation 
given by Prugniel & Simie^ (1 19971) 

j(r) = ji](r/a s ), 

j(x) * jc^exp(-x 1/m ), 

. = f F(2m) } Ip 

Jl \r[(3-p)m]/2fl; 

p a 1 -0.6097 /m + 0.05463 /m 2 , 

where the latter equation comes from iLima Neto et a fl (fl999h . 
F(x) is the Gamma function, and Iq, a s , an d m refer to the 
Sersic p arameters from Eq.[2] We also refer to lMamon & Lokasl 
(2005b) for a compilation of the relevant formulae. 

Integrating the luminosity density, j(r), the enclosed lumi- 
nosity is found. ILima Neto et alJ (jj.999) give 
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Fig. 7. A comparsion between the enclosed luminosities of the 
double /3 profile (dashed blue line) and the double Sersic profile 
(red solid line), normalized to the double Sersic total luminosity. 
The arrow indicates the last spectroscopic point measured on the 
major axis. 



L s (r) = L tot L s (r/a s ), 

y[(3-p)m,x'/ m ] 
LAx) — , 

r[(3 -/>>»] 

where y is now the incomplete Gamma function and 
Ltot = 2mnT(2m)I()a 2 



(7) 
(8) 



(9) 



is the total luminosity of the profile. For our double Sersic pro- 
file, the luminosity is simply the sum of the luminosity of both 
the inner and outer Sersic profiles. 

Since beyond ~55 kpc the double-/? and Sersic profiles dis- 
play the most distinct behaviors (cf. Fig. [3J, we inspected the 
influence of these dissimilar extrapolations on their respective 
enclosed luminosities, which are directly related to the stellar 
mass through the mass-to-light ratio. Fig. Q considers the en- 
closed luminosities of both models normalized to the total lumi- 
nosity of the double Sersic model given by Eq. [9] As expected 
from Fig. [3] the luminosities are very similar not only within the 
~55 kpc where the galaxy light was directly measured, but well 
beyond. A ~10% difference in the total enclosed luminosity is 
reached only at 200 kpc. In practice, however, the extension of 
elliptical galaxies may be lower than this limit. In NGC 4636, 
for exam ple, the globu lar cluster distribution is cut-off at about 
-60 kpc dPirsch et alj |2005). A similar behavior has been seen 
in M 87, the central galaxy of Virgo, where the planetary nebu- 
lae distribution show a truncation of the stellar halo at ~150 kpc 
dPohertv et alj |2009). Therefore, any differences in the extrapo- 
lations beyond 200 kpc are not relevant to a dynamical analysis 
of the galaxy. 

Given the similarities between the light profiles and to avoid 
any uncertainties introduced by the approximations used in the 
Sersic profile deprojection (which are not present in the analytic 
form of the /3 profile deprojection), the Jeans analysis that we 
now describe was based on the stellar distribution given by the /3 
model profile. In any case, the velocity dispersions shown in Fig. 
[8] were found to vary by less than 2 km s at the most distant 
kinematic point, when the double Sersic model was used a basis 
of the stellar mass profile. 



4.2. Jeans analysis 

The basis of our modeling is the spherical Jeans equation, which 
is strictly valid only for infinite, non-rotating systems. Its so- 
lutions, when applied to real galaxies, may be considered as 
good approximations, but there is no guarantee that they cor- 
respond to physically correct distribution functions. Moreover, 
the potential and the orbital an i sotrop y are degenerate, although 
we refer to lHansen & Moord (|2006) and the remarks below. 
On the positive side, this equation has an appealing simplicity, 
known anisotropies of elliptical galaxies in the inner regions are 
generally small (e.g. Ger hard et al.ll200 ll ICappellari et al . 2006) 
and one can use the higher mome nts of the velocity dis tribu- 
tion as constraints of the anisotropy. Cappell ari et al.l (2006) also 
demonstrated that more sophisticated dynamical modeling does 
not reveal grossly different results. The main uncertainty proba- 
bly does not arise from the modeling approach, but collectively 
from the distance uncertainty, the quality of the data, and the 
assumption of sphericity. 

For convenience and clarity, we give in the following the 
basic formulae. The Jeans equation is (e.g. Binnev & Tremaine 
2008) 



d[j(r)<r 2 r (r)] 2/3(r) , 
j + l(r)(K 



r r l 



(10) 



which relates the light distribution, j(r), and radial velocity dis- 
persion, cr,.(r), with the underlying gravitational potential or en- 
closed total mass, M(r), and where 



(3(r) = 1 - o-l/o-l 



(ID 



is the anisotropy parameter, which indicates possible departures 
from pure isotropy in the form of radial (J3 > 0) or tangential 
08 <0) orbits. 

For a given light and mass distribution, and a constant 
anisotropy, the solution to the Jeans equation is 



j(r)(r, 



(12) 
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(e.g. |M 

amon & Lokas 2005b, hereafter ML05), which is then 
projected to be compared with the velocity dispersion measure- 
ments 



LOS 



(R) = 



2 

W) 



r r 

Jr Vr 2 - R 2 Jr 



Pjv r dr 



R 2 



(13) 



dBinnev & MamonTl 982). For the anisotropy profiles adopted in 
Sects. l4~4l and 14.51 Eqs. [T2land[T3lcan be combined to yield the 
single integral solution 



2 r°° dr 
"los(*> = T^j J r K(R,r)j(r)M(r)-, 



(14) 



where the kernels K for t he different anisotropy m odels are pro- 
vided in Appendix A2 of (Mamon & Lokas 2005b). 

The accuracy of the models can be quantified by the use of a 
simple x 2 test 



X 



-model \2 



6cr° 



.obs 



(15) 



where o" obs is the measured velocity dispersion, cr model is the pre- 
dicted dispersion at the same radial distance, 5cr obs is the obser- 
vational error associated with the dispersion measurement, and 
the sum is taken over all radial bins. 

Under the spherical Jeans analysis, additional constraints 
on the dynamical behavior of the galaxy can be obtained by 
modeling the higher order moment of the velocity distribu- 
tion which can be related to the Gauss-Hermite coefficient /i 4 
(Ivan der Marel & Franx| [l993). For this, the higher order Jeans 
equations must be solved. This task can be achieved only if ad- 
ditional information about the distribution function is known. 
In particular, by assuming distribution functions of the form 
f(E, L) = f (E)L~ 2 P with constant /3, it is possible to reduce 
the two higher-order Jeans equation to only one. A more de- 
ta iled discussio n, as well as the relevant for mulae can be seen 
in lLokasl(l2002l) and iNapolitano et all (120091) . Our results using 
this approach are presented and discussed in Sect. 14.61 

4.3. The velocity dispersion profile 

To construct the radial dependence of the velocity dispersion 
profile, we averaged for a given radial bin the dispersion along 
the major and the minor axis (although the ellipticity is tiny) 
from Figs.[5]and[6] The dispersion values are symmetrical with 
respect to the galaxy center and are in agreement between the 
both axes when errors are considered. The bins are chosen in 
such a way that the inner region, where the S/N is highest, has 
the densest binning while the outer radii are covered with larger 
bins. The principal uncertainties are those in the mean values cal- 
culated by pPXF Final values are listed in Table [4] which also 
lists the adopted /14 values. Their uncertainties are the standard 
deviations of the /14 values in the respective bins. 

4.4. Models without DM 

The simplest model is the one in which M(r) = M*(r) (i.e. no 
DM) and /3 — (orbital isotropy). One can adopt the expected 
value for the M/L from an stellar populations analysis (which 
is a hypothetical value) or treat it as a free parameter in the 
Jeans equation solution, which is the procedure we chose here. 
This gives an upper limit to the stellar M/L. This model is 
shown as the black solid line in Fig. [8] (top left panel). It uses 



Table 4. Adopted velocity dispersion,cr, and /14 values and their 
respective uncertainties. We refrain from giving the ^-values 
since they are not modelled. 



Radius 


a 


Ao- 


hn 


Ah 4 


(arcsec) 


(km s _1 ) 


(km s" 1 ) 






0.5 


239 


2 


0.051 


0.013 


1.5 


232 


2 


0.045 


0.013 


2.5 


222 


2 


0.010 


0.017 


3.5 


217 


2 


0.005 


0.018 


4.5 


211 


2 


0.002 


0.011 


5.5 


205 


2 


0.055 


0.014 


6.5 


202 


2 


0.043 


0.027 


7.5 


199 


2 


0.042 


0.016 


8.5 


189 


2 


0.043 


0.019 


9.5 


193 


2 


0.033 


0.012 


10 


191 


2 


0.040 


0.015 


13 


190 


2 


0.048 


0.026 


16 


181 


2 


0.037 


0.026 


22 


167 


2 


0.019 


0.023 


28 


153 


2 


0.042 


0.036 


42 


141 


3.3 


0.044 


0.033 


52 


135 


6.7 


0.008 


0.022 


62 


131 


14.3 


-0.012 


0.036 


85 


125 


18.7 


-0.013 


0.015 



an M/Lr- 3.13 and already gives a very good description of 
the dispersion profile with^ 2 =46.8 for 16 degrees of freedom, 
where the statistics is calculated outside 400 pc, which is the 
radius where the data are no longer affected by seeing. 

When the anisotropy was allowed to vary, we obtained al- 
most the same M/Lr of 3.11 and /? = 0.12, which is a slightly 
better fit (red dotted line). Table [5] shows the parameters for dif- 
ferent best-fit models. An unrealistic totally radial model (J3 = 1) 
is also shown in Fig.[8]as a comparison ( blue dashed line). 

From theoretical considerations (e.g. Abadi et al. 20061), el- 
liptical galaxies are expected to have an anisotropy profile that is 
increasingly radially biased with growing radius. An analytical 
profile with this characteristic was introduced by ML05 as 



PmlV) = fo- 



ri- r„ 



(16) 



with fio ~ 0.5 and r a ~ 1 A R e when fitted to the merger simula- 
tions of lDekel et alJ (120051) . Another pr ofile frequently used is 
the O sipkov-Merritt anisotropy model (Osipkovl ll979t iMerrittl 



the O sipk 
119851) . 

Poiu(r) 



(17) 



where r a is a characteristic radius marking the transition to outer 
radially biased orbits in both cases. Since the data do not re- 
ally constrain the value fo r r a , w e fixed it to r a - l.4R e , using 
R e -3l" from Faber et al. (1989). Even though our derived ef- 
fective radius is significantly larger (see Sect. 13. U . we use this 
older estimation since it has been used in the previous studies 
of the galaxy, and also because it implies a faster transition to 
a radial behavior. The predictions from the Jeans equation us- 
ing these anisotropies are shown in Fig. [8] (top left panel). The 
ML05 anisotropy gives a good description for the dispersion pro- 
file (green dashed line), but formally not better than the con- 
stant /3 = 0.12 model. It reaches a value of (3 — 0.33 at the outer- 
most kinematic point. However, the Osipkov-Merritt anisotropy 
(green dash-dotted line) overpredicts the dispersion for the in- 
ner regions and underpredicts it for larger radii, making it less 
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Fig. 8. Dynamical models of NGC 7507. Black circles indicate in all panels the line-of-sight velocity dispersion. All models are 
given outside of 400 pc, where the luminosity model is unaffected by seeing. Top left: Stars-only models. The black solid line is an 
isotropic model using M/Lr = 3.13. This model is shown in all panels as a comparison. The red dotted line is the best-fit model for 
a constant anisotropy using f3 = 0. 12, while the blue dashed line is a fully radial model (J3 = 1). The green dashed line is the model 
using the anisotropy as given in Eq.[l6] while the gr een dot-dashe d line is a model using the Osipkov-Merrit anisiotropy (Eq.fTTIi. 
The red solid line is a model with M/Lr-2, based on lTortora et alj 's stellar population analysis (see Section lBTTl i. Top right: MOND 
predictions using the "simple" and "standard" interpolation formulae. The dashed red area is an isotropic model, while the blue area 
indicates a fully radial model for the former interpolation. Bottom left: NFW models using c — 5 and c — 10. Parameters for each 
model are reported in Table [6] Bottom right: Logarithmic dark halo models. Using /3- 0, f3 — 0.3, and ML anisotropy. Parameters 
for each model are given in Table [7] 



Table 5. Comparison of different best-fit dynamical models 
without DM. Errors indicate the 1 - <x confidence levels. 



Model 




P 


x 1 




W Q /L QM ] 






Stars only 


Q 1O+0.01 

-Q-Q2 


0.0 


46.8 


Stars only + est. p 


3 i j +0.02 


12 +008 


44.1 


Stars only + ML" 


3 01+^ 


0.33 


54.5 


Stars only + HM 6 


9 qo+o.02 

y -0.02 


0.45 


90.5 



Notes. (fl) Indicates the ML05 anisotropy. The value for /? is the max- 
iumum value for the anisotropy reach ed at the last kinematic point. 
ib> Indicates the Hansen & Moore ( 2006) anisotropy modelled using Eq. 
[I6]with r a = 500 pc. The value of /? has the same meaning as before. 



preferable as a realistic description within the context of this 
model. 



An interesting relation between anisotropy j3 and the loga- 
rithmic density slope a of a spheroida l system near equi l ibrium 
has emerged from DM simulations. Hans en & Moord ([2006) 
found /3(a) = 1-1.15(1+ a/ 6), a lbeit with some scatt er. This 
relation has also been analyzed by M amon et alJ (120061) for el- 
lipticals formed in simulations of mergers of spirals consisting 
of stars, gas and DM. They found different results for the DM, 
but a quite similar relation for the stellar density distribution. If 
this were universal, the dynamical analysis of spherical systems 
would be greatly facilitated. We apply the above relation to NGC 
7507 and calculate for our double-beta model (Eq. |4]l the loga- 
rithmic slope a = (r/j)dj(r)/dr. The resulting anisotropy pro- 
file shows a swift increase from isotropy at the center to /3 = 0.3 
at around 500 pc, with a mild increase outwards. This behav- 
ior can be reasonably reproduced using the ML05 anisotropy 
profile with a small r a . In the context of a no-DM model, this 
anisotropy favors an even lower stellar M/L and underestimates 
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the observed velocity dispersion at large radii, giving a poor fit 
to the data. 



4.5. Including a dark halo 

Even though the velocity dispersion profile up to ~ 90" of 
NGC 7507 can satisfactorily be described by models without 
dark matter, the paradigm of galaxy formation requires the pres- 
ence of dark halos. We there fore allow for a (cuspy) NFW pro- 
file (e.g. Nava rro et al.l lT996) and a (cored) logarithmic halo. In 
both cases we fix the stellar M/L to 3.01, therefore maximiz- 
ing the contribution of the stellar component to the tot al mass 
budget ("maximum spheroid", e.g. lWeiimans et aD2009l) . this is 
supported by the fact that the contribution of D M in the inner 
regions is thought to be close to negligible (e.g. | Gerhard et all 
2001; lBorriello et alJl2003l:lMamon & L okas 2005a|). 



4.5.1. NFW models 



The NFW profile is given by 

Ps 



PnfwM 



(r/r s )(l + r/r s f 



(18) 



where p s and r s are a characteristic density and radius, respec- 
tively; the enclosed mass for this profile is given by 



M NFW (r) = 4 W ] I In ( 1 + r/r s ) - ' ' ' 



1 + r/r s 



(19) 



We explored three models, an isotropic model, a model with 
a slight radial anisotropy (y8=0.3) and a model using the ML 
anisotropy with r a - 43" as in Sect. 14.41 We excluded the 
possibility of tangential orbits beacuse they were inconsistent 
with the stars-only models. For each anisotropy, we used a set 
of concentrations, c = {1, 1.25, 1.5... 13), which encompasses 
the rang e of concentrations fou nd for low luminosity ellipticals 
(c ~ 5. iNapolitano et aLll2008l) and the on es found in simula- 
tions for halos of similar mass (c ~ 10, e.g. lBullocketal.ll200lt 
iMaccio et al . 2008). Concentrations are related to the character- 
istic density by 



200 



Pcrit 



3 ln(l +c)-c/(l +c)' 



(20) 



where the critical density is p CI j t = 3H^/%nG. We also set a grid 
of characteristic radii, r s = {1, 2, 100) kpc, and calculate 
the x 2 f° r eac h ( c > r s) P a ir- The results can be seen in Fig. [9] 
68%, 90% and 9 9% confide nce levels correspond to Ay 2 = 
{2.30, 4.61, 9.21} <|Ayjij|ll976|). 

As seen in Fig.|9j we cannot constrain simultaneously c and 
r s , with the available data/modeling. Even though we have ex- 
plored a generous range of parameters, a global minimum lies 
beyond the explored space. The minimum x 1 values for the 
isotropic and ML models are worse than the best-fit stars-only 
models, even though the models with an NFW profile include 
more free parameters. This is an indication that the data shows 
a preference for the stars-only models, or more exactly, that our 
assumed stellar M/Lr is closer to the total M/L, which should 
be constant across the galaxy. The exception is the model with 
y6 = 0.3, which gives a significantly smaller^ 2 than all of the no- 
DM models. 

Since we cannot constrain simultaneously the values for c 
and r s , we explore the amount of DM in the NFW halos by se- 
lecting the cases c — 5 and c — 10. The NFW parameters of each 



halo and the assumed concentrations can be seen in Table[6] The 
velocity dispersion prediction for each of these models is de- 
picted in Fig. [8] (bottom left panel). 

The best-fit model is the one for c-5 and = 0.3, which 
gives a virial mass of 2.2 x 10 11 Mq. The c=10 model for 
the same anisotropy gives a slightly worse fit and a signifi- 
cantly lower virial mass with 5.9 x 10 10 Mq. Even though the 
ML anisotropy should give a more realistic representation of 
the true orbital structure of elliptical galaxies, the mass which 
correspond to the best-fit scale radius is very low, although the 
1 — <t level i s broad with a maximum mass of 1.9 x 10 10 Mq 
for c — 5. The lHansen & Moord (l2006) anisotropy performs even 
more poorly, with even lower predicted masses and much larger 

x 1 - 

The fraction of DM within \R e (using the value from 
Faber et al . 1989) is lower than 5% for all models. This number 
can be compared wit h the values found for galaxies of similar 
brightness studied by IWeiimans etaf] (120091) : NGC 3379 with 
fDM(<Re) = 0.08 and NGC 821 with f DM (<R e ) = 0.18. Even 
though this low value could be considered as the natural result of 
the minimal halo assumption, this is not exactly the case. If we 
were to drop this assumption and allow the stellar mass-to-light 
ratio to be a free parameter for the c — 5 and c — 10 NFW ha- 
los, we would find that the M/L diminishes at most to a value of 
2.9, depending on the assumed anisotropy, reinforcing the prior 
assumption. DM fractions within 2R e can also be seen in Table 
i 



4.5.2. Models with a logarithmic halo 

Even though most, if not all, N-body simulations favor cuspy 
profiles, observations in dwarf and low surface brightness galax- 
ies indicate the presence o f constant density cores instead of 
r 'cusps (see lde Blokl2010l for a recent review). From the differ- 
ent types of cored profiles we selected the logarithmic potential 
(Binn ev & Tremainell2008l) . which has a density 



PlogW = 



3r 2 + r 2 



47rG(r 2 + r 2 ) 2 ' 
and enclosed mass 



M h0G (r) = 



v 2 r 



(21) 



(22) 



The modeling procedure is similar to the one applied in the 
previous Section. As for the NFW profiles, the stellar M/Lr 
remains fixed at 3.01. A grid of parameters was then set with 
r = {1,2. ..80} kpc and v = {40, 45. ..215} km s~' and for each 
pair, the value ofx 2 is calculated using Eq.[l5]for the same three 
anisotropics as in the previous section; /3 = 0, /? = 0.3, and j3(r) 
following Eq. [16] The grid was later refined for the case /3 = 0. 
The best-fitting halo parameters for each studied case are given 
in Tabled 

The isotropic case allows a light halo with a very small char- 
acteristic radius of only 25 pc. Assuming/3 = 0.3 the halo is much 
larger with ro = 7 kpc and a circular velocity of 100 km s . This 
radial model provides a significantly tighter fit to the observed 
velocity dispersion. This halo is the one that gives the lowest^ 2 
amongst all models (with or without DM). The ML anisotropy 
allows a large range of possible halos, but giving a worse fit than 
the p = 0.3 case. As in the case of NFW halos, the amount of DM 
enclosed within one R e is minimal, again having values of less 
than 5% for all models. The velocity dispersion for these three 
best-fit models can be seen in Fig. [8] (bottom right panel). 
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Table 6. Comparison of different best fitting dynamical models including NFW halos. 



P ps r s ^200 M(r 2Kl )IL R /dm < R e /dm < 2R e 

[Mqpc- 3 ] [kpc] [kpc] [M ] [M /L ] [%] [%] 

Stars (Y R =3.01) + NFW (c=5) 

0.0 1.191 x 10 3 30 150 3.87+" x 10" 7.3 3.4+"| 7.5+" 63.0 

0.3 1.191 x 10" 3 25 125 2.24« : ° x 10" 5.6 2.7+° : jj 8.9+f| 40.5 

ML 1.191 x IP 3 1 5 1.4 +ll? xl0 7 3A <1_ < j 54.5 

Stars (Y Ji =3.01) + NFW (c=10) 

llO 6.135 x 10- 3 9 90 8.36+^ x W> 4~1 JJ^l TA^ 593 

0.3 6.135 x 10- 3 8 80 5.87+g x 10 10 3.8 3.3;*! 8.6+ 22 42.9 

ML 6.135 x 10 3 1 10 0.ir"xl0 9 3.1 < l' <1 54.5 



Notes. Column 1 indicates the value for the anisotropy considered for the stellar population. ML indicates the Mamon & Lokas ( 2005b) anisotropy 
profile with parameters decribed in the text. Column 2 gives the characteristic density of the NFW halo, which comes from the chosen concentration 
as given in Eq. [20] Column 3 gives the best fitting scale radius for the chosen concentration. Columns 4 and 5 give the virial radius and mass of 
the DM halo, respectively. Masses are also indicated with their 1-cr confidence l evel. Column 6 i ndicates the virial mass to light ratio. Columns 7 
and 8 give the DM fraction within 1 and 2 R e (using the effective radius given by Faber et al. 1989), defined as / DM = M NFW /(M NFW + M t ), where 
the baryonic mass is given by Eq.[6] Column 9 gives the^ 2 for each model. 
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Fig. 9. x 1 for models including an NFW halo. Minimum x 1 f° r eacn model is indicated. Left panel: NFW halo and stars with 
isotropic distribution. Central panel: NFW halo plus stars with a radial anisotropy of /3-0.3. Right panel: NFW halo and stars 
following a ML anisotropy profile. 



Table 7. Comparison of different best-fit dynamical models us- 
ing a logarithmic DM halo. All models use M/L* « = 3. 



p 


ro 


Vo 


fDM < Re 


/dm < 2/? c 


x 1 




[pc] 


km s -1 


[%] 


[%] 




0.0 


25 


60 


4.1 


5.4 


47.7 


0.3 


7 000 


100 


2.3 


7.4 


38.6 


ML 


74 000 


210 


< 1 


< 1 


52.4 



tions since they would lie between the two aforementioned mod- 
els, without adding new information. As in the cases of the ve- 
locity dispersion profiles, since the light profile is affected by 
seeing, the model is not expected to reproduce the inner 400 pc. 



4.6. The Gauss-Hermite parameter h A 

Predictions for the Gauss-Hermite /14 profile for different models 
can be seen in Fig.[l0l where we have followed the procedure of 
iLokasI (12001 and related kurtosis and h\ using 

k(R) * 8 V6/24, (23) 

(I van der Marel & Franxl [19931) . We have considered the stars- 
only, NFW with c = 5 and 10, and the logarithmic cases 
withr isotropy and completely radial anisotropies. Although 
the velocity dispersion profiles are best reproduced with radial 
anisotropies using /3 = 0.3, we did not calculate their /14 predic- 



The measured /14 data are noisy (green circles) and all mod- 
els system atically fail to reproduce its high values. We note that 
the iK+OOl data are equally noisy and reach similar high values 
(blue stars in Fig. fTOb. The d ifference is most evident in the in- 
ner 500 pc where the IK+OOl values are closer to zero and ours 
scatter around ~0.05. Although some systematic effect may be 
shifting our measurements to higher values, it can hardly be S /N 
(which is especially high at small radii) nor template mismatch- 
ing (innermost spectra are very well-fitted using the comprehen- 
sive MILES stellar library). Perhaps the only conclusion that we 
can give is that the ^4 data does not favor isotropic models, which 
predict slightly negative values for /14 (see Sect. I5.3l for an alter- 
native explanation). 
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Fig. 10. Kurtosis models for NGC 7507. Green circles repre- 
sent our hi, measurements as given in Table [4] Blue stars are 
the measurements from lK+OOl Black, red, blue, and green lines 
represent the stars-only, NFW (c=5), NFW (c= 10), and logarith- 
mic models, with parameters described in the text, respectively. 
Solid lines are isotropic models, while dotted lines are fully ra- 
dial models. 



5. Discussion 

5.1. Stellar M/L 

The stellar M/L determines the amount of baryonic matter and 
thus, based on its difference from the total mass, also the fraction 
of dark matter. Our M/L-value is not quite consistent with what 
one would theoretically expect for an old, metal-rich population 
by comparison with published population synthesis models, but 
it i s in agreement with literature values. 

IK+OOI gives a central M/Lb value of 5.9 for a distance 
of 26.9 Mpc. Since dynamically derived M/Ls are inversely 
proportional to distance, at our adopted distance of 23.2 Mpc 
this value transforms i nto M/Lb = 6.8, wh ich gives M/Lr = 3.14 
using (B - R) = 1 .90 (iFranx et al.1 [l989bl) and (B-R)q = 1.06 
(Binnev & Merrifield 1998), in excellent agreement with our 
value. Magorrian & Ballantvne (2001), based on kinematic data 
bv lBertin et alJ (1 19941) . found M/Lg ~ 7, somewhat higher than 
our value, but also in agreement considering their shorter dis- 
tance (21 Mpc). 

The stellar M/L s even of old ell i ptical galaxies are appar- 
ently not universal: Cappell ari et al ] hood) find values in the 
range 1.5 < M/L/ < 6.0. For NGC 7507, M/L s = 3.1 trans- 
lat es into M/Lj - 2.42. If w e only consider elliptical galaxies in 
the ICappellari et al.1 d2006l) sample, that is, excluding SO galax- 
ies, this range would narrow to 2.3< M/L? <6.1, putting NGC 
7507 among the ellipticals with the lowest M/L. This value, 
unless the metal abundance of NGC 7507 is distinctly subso- 
lar, would indicate a somewhat younger galaxy, with an age 
around 8 Gyr for a solar abundance, which a grees with previous 
line-strength measu rements of this galaxy dTrager et al.l [l998: 
lOgando et al. I I20081) . This is unsurprising because a large frac- 
tion of isolated elliptical galaxies are kn own to be of younger age 
(ICollobert et alJ2006l:lReda et al.l2007|). which is c onsistent also 
with simulations (iNiemi et al.ll2010l) . lTortora et all 's (2009) stel- 
lar population models give a very similar age of 7.7 Gyr (Kroupa 



IMF), but a siginificantly different M/Lr = 2. This value is in 
contradiction with the line strength analyses and with our dy- 
namically derived M/L and would imply a huge DM content, 
even in the innermost parts of the galaxy (red solid line in Fig. 
[8] top left panel). 



5.2. Dark matter? 

That an isotropic model with no DM closely fits the observa- 
tions nicely, does not imply that there is no DM. This situa- 
tion was realized before in analyses where PNe were used as 
dynamical tracers (Romano w skv et all 120031: iDekel et all [2005. 
ML05). However, the relatively low stellar M/L-value of NGC 
7507 is in line with the finding that, regardless of the situation 
at larger radii, the DM content in the inner region of e llipti- 
cal galaxies is negligible (|Gerhard et al.ll2.00ll: iMamon & Lokasl 
2005a; Cappellari et al. 2006). Within the framework of a spher- 
ical model, a possibility to host DM at larger radii is a radial 
anisotropy. It must then be a coincidence that the observed dis- 
persion so closely follows the predictions of the photometric 
model under isotropy. This is a difference to the use of PNe as 
dynamical tracers, where the density profile of the parent popu- 
lation is notoriously uncertain. 

Are the stellar orbits of NGC 7507 affected by a strong radial 
anisotropy? The comparison with the kurtosis profiles for a few 
models shows that a fully radial anisotropy of the stellar compo- 
nent can produce a consistently positive kurtosis of the observed 
order, although failing in the innermost regions. However, we 
say that wit h caution, since the kurtosis might also be influenced 
by rotation dDekel et al.ll2005l) . 

NGC 7507 with a stellar mass of about 2 x 10 11 M Q belongs 
to the most ma s sive i solated elliptical galaxies in the simulations 
of Niemi et al.l (l2010h . The dark halos in which they are embed- 
ded have predicted masses of the order of a few times 10 12 Mq 
within their virial radii. Our halos from Table|6]are far from this 
mass, although we cauti on that the extrapolation to the virial ra- 
dius might be uncertain (M amon & Lokas 2005^). 



5.3. Limitations of a spherical model 

Even though the apparent spherical symmetry of the galaxy and 
the small amount of rotation give reasons to use our modeling 
approach, at the same time other observations may point to a 
more complex scenario. High /14 values could be signature of 
an undetected face-on disk or flattening along th e line of sight 
( Mag orriarJl 1 9991: iMagorrian & Ballantvne 2001), which might 
also explain t he failure to reproduce the h\ values when using the 
lLokasI d2002l) simplification. The second observed particularity 
is the rotation along the mino r axis. E lliptical galaxies with this 
property are rare: IFranx et all dl989bl) found that minor axis ro- 
tation was larger than major axis rotation in only 2 out of 22 
galaxies in the ir sample. Th is behavior is a signature of triaxial 
galaxies (e.g. lBinnevlll985l) . Another sign of triaxiality comes 
from the outermost velocity along the major axis which shows a 
significant difference with the inner measurements. This change 
in the mean velocity corresponds to a sudden change in the el- 
lipticity at the same radius (~100", Figs.|2]and|5j. 



5.4. MOND and the baryonic Tully-Fisher relation 

Among the problems that MOND faces is how the Newtonian 
and MONDian regimes are linked. Taking this interpola- 
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tion function as JJ-(x) — x/(l + x) (known as "simple", 
Famaey & Binnev 2005), the MONDian circular velocity is 



V 2 „ I U 4 XT 

2 . . circ.N / circ,N 9 ,^*\ 

"circMONDM = ~Y~ + \ + a ° U ci^N r ( 24 ) 

(e.g. iRichtler et al.l 120081) . where 0ctrc,N is the Newtonian cir- 
cular velocity associate d with the stellar co mponent and ao = 
1.35+jJ^ x 10 _8 cms~ 2 dFamaev et al.1 120071) is the acceleration 
constant that separates the MOND and Newtonian realms. An 
alternative interpolation function, dubbed "standard", is /z(x) = 
x/ Vl +x 2 dSanders & McGaughl 120021) . for which the circular 
velocity is given by 



v 4 v & 

4 , , circ,N circ,N 9 4. o 

y circ,MOND('-) = + V + fl U circ,N^ < 25 ) 

(e.g. Samurovic & Cirkovic 2008). The masses associated with 
these circular velocities are introduced in Eq. [14] to obtain the 
MOND predictions for the projected velocity dispersion profile 
for the isotropic case, using in all cases M/L+ r = 3. As seen in 
Fig. [8] (top right panel), the simple interpolation formula cannot 
explain the velocity dispersion profile under isotropy. The out- 
ermost measured points can only be explained under a totally 
radial model and only when the lowest allowed values for ao 
are considered. The standard interpolation formula (green lines 
in the same plot) gives very similar predictions for the isotropic 
and radial cases for the outermost bins, and again favors a lower 
value for ao. 

The barred spiral galaxy NGC 7513 is found at an angular 
distance of 18.2 arcmin from NGC 7507. In the MOND context, 
it is interesting to ask whether it could provide an external grav- 
itational field which would modify the MOND phenomenology 
in NGC 7507. The mean of the published Tully-Fisher distances 
is 19.3 Mpc and thus places NGC 7513 with respect to NGC 
7507 (23 Mpc) somewhat in the foreground and a separation of 4 
Mpc would produce a very weak field. Adopting an inclination- 
corrected Hl-line width of 275 km s" 1 as the valu e for a constant 
MONDian rotation curve (iSpringob et al.1 12007). the accelera- 
tion at NGC 7507 would be 6.3 X 10~ 13 m s~ 2 . Assuming the 
same distance as for NGC 7507, the separation would be 123 
kpc and the acceleration at NGC 7507 would be 2 x 10~" m 
s -2 . Therefore, an external field should not influence any of the 
MONDian dynamics of NGC 7507. 

The baryonic Tully-Fisher relation (BTFR) of spiral galax- 
ies relates circular velocity (which is a proxy for the total mass) 
to the baryonic content of the galaxy. This relation, which reads 
Mbar = 50 V 4 , where uj is the circular velocity in the flat part of 
the rotation cur ve, has been found to be correct over ma ny orders 
of magnitude (lMcGaughl2005UTrachternach et al.l2009l) . A tight 
relation between baryonic mass and circular velocities is not ex- 
pected when DM is assumed to dominate t he dynamics, but it 
finds a natural explanation under MOND dSanders & McGaugh 
120021) . Similarly to spirals, early-type galaxiess have been found 
to follow a relation that is analogous to the Tully-Fish er relation 
dGerhardet alJ '2001: Magorrian & Ba llantvndl2001l) . but since 
circular velocities up to a large radius are difficult to obtain, there 
is no surprise in finding that ellipticals follow a relation that is 
offset from the spiral BTFR, i.e. circular velocities are probably 
overestimated. For NGC 7507, we calculated the circular veloc- 
ity at the outermost point with measured kinematics of 250 km 
s -1 . Using McGaugh's version of the BTFR, this implies a mass 
of 2xlO n M©. Using Eq.|6]with our preferred Mj L value, we 



obtained 2.5xlO n M , that is, NGC 7507 would lie somewhat 
above the BTFR for spirals; this appears inconsisten t with pre- 
vious results for other ellipticals dGerhard et alJ uOOl). although 
consistent when all uncertainties are taken into account. We also 
note that the sl ope of 4 required b y MOND remains a disputed 
value (see e.g. lGurovich et al.ll2010l) . 

6. Summary and conclusions 

We have obtained wide-field photometry in Kron-Cousins R and 
Washington C as well as new long-slit spectroscopy of the field 
elliptical NGC 7507 out to about ~ 90", reaching farther out 
than previous studies. 

We have measured almost no rotation along the major axis 
and significant rotation (Av ~ 50 kms -1 ) along the minor axis 
(although the galaxy is almost perfectly circular). The velocity 
dispersion profile of the galaxy shows a rapid decline along both 
axes. 

We have performed a spherical Jeans analysis to find the 
mass profile that most closely represents the projected veloc- 
ity dispersion profile. When assuming isotropy, a radially con- 
stant M/Lr ratio of 3.1 would be sufficient to reproduce the 
velocity dispersions. This value corresponds, if solar abun- 
dance is assumed, to an age in the range 8-10 Gyrs. DM ha- 
los can be included in the isotropic case, but provide signif- 
icantly poorer fits than the stars-only models. When we al- 
lowed for radial anisotropies and cored/cuspy dark halos, we 
were able to improve the fit to the velocit y dispersion pro- 
file ma rgin ally. This contradicts the result od Krona witter et al.l 
(120001) and iGerhard et ail d200ll) who found this galaxy to be 
one of the best examp les for hosting a dark halo . It is h owever in 
line with the results of lMagorrian & Ballantvne (2001) who also 
found no clear evidence of DM. 

The most massive NFW halo, using j3 = 0.3 and a concen- 
tration parameter of 5, has a scale r s = 25 kpc, implying a virial 
mass of only 2.2 +2 y l 2 x 10" Mq. Modified Newtonian dynamics 
(MOND), if applied straightforwardly, predicts velocity disper- 
sions that are too high. 

It appears that NGC 7507 is a very interesting case in the 
discussion of DM in elliptical galaxies. It remains to be seen 
whether non-spherical models based on an extended data set, 
perhaps including also PNe and/or globular clusters, would pro- 
vide different results. 
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Table 1. T\ CTIO/MOSAIC surface photometry of NGC 7507. First column indicates the semi-major axis from the ellipse fitting. 



R 




6 


fl 4 


PA 


(arcsec) 


(mag arcsec -2 ) 









0.148 


15 


.136 


± 0.001 


0.264 


0.104 


± 0.043 


-23.9 


± 8 


.4 


0.511 


15 


.186 


± 0.001 


0.052 


-0.023 


±0.011 


-33.6 ± 16 


.3 


1.095 


15 


.550 


± 0.001 


0.031 


-0.003 


± 0.001 


8.3 


± 3 


.4 


1.457 


15 


.797 


± 0.001 


0.033 


-0.001 


± 0.001 


8.8 


± 2 


.1 


2.133 


16 


.194 


± 0.001 


0.040 


-0.001 


± 0.001 


13.8 


± 1 


.5 


2.581 


16 


.424 


± 0.001 


0.043 


-0.001 


± 0.001 


12.8 


± 1 


.3 


3.123 


16 


.678 


± 0.001 


0.040 


-0.002 


± 0.001 


15.4 


± 1 


.0 


3.435 


16 


.810 


± 0.001 


0.039 


-0.001 


± 0.001 


16.3 


± 1 


.2 


4.572 


17 


.200 


± 0.001 


0.046 


-0.000 


± 0.001 


14.4 


± 1 


.1 


5.029 


17 


.326 


± 0.001 


0.048 


-0.000 


± 0.001 


16.4 


± 


.8 


5.532 


17 


.449 


± 0.002 


0.052 


-0.001 


± 0.001 


15.1 


± 1 


.1 


6.086 


17 


.574 


± 0.002 


0.054 


-0.002 


± 0.001 


15.7 


± 


.9 


7.364 


17 


.837 


± 0.002 


0.054 


0.003 


± 0.001 


15.6 


± 


.8 


8.100 


17 


.983 


± 0.003 


0.050 


0.001 


± 0.001 


16.7 


± 


.8 


8.910 


18 


.138 


± 0.003 


0.046 


0.001 


± 0.001 


14.7 


± 


.9 


9.801 


18 


.307 


± 0.004 


0.042 


-0.001 


± 0.001 


15.5 


± 


.9 


10.782 


18 


.475 


± 0.001 


0.042 


0.002 


± 0.001 


12.1 


± 


.5 


11.860 


18 


.646 


± 0.001 


0.042 


0.002 


± 0.001 


12.1 


±0 


.5 


13.046 


18 


.817 


± 0.001 


0.042 


0.003 


± 0.001 


12.1 


± 


.5 


14.351 


18 


.985 


± 0.001 


0.042 


0.003 


± 0.001 


12.1 


± 


.6 


15.786 


19 


.145 


± 0.001 


0.045 


0.002 


± 0.001 


13.0 


± 


.5 


17.364 


19 


.302 


± 0.001 


0.048 


0.001 


± 0.001 


11.8 


± 


.5 


19.101 


19 


468 


± 0.002 


0.043 


0.004 


± 0.001 


10.3 


± 


.6 


21.011 


19 


.621 


± 0.002 


0.045 


0.003 


± 0.001 


10.3 


± 


.5 


23.112 


19 


.772 


± 0.002 


0.047 


0.003 


± 0.001 


8.4 


± 


.5 


25.423 


19 


.928 


± 0.003 


0.046 


0.002 


± 0.001 


9.2 


±0 


.5 


27.965 


20 


.096 


± 0.003 


0.044 


0.001 


± 0.001 


10.3 


± 


,4 


30.762 


20 


.278 


± 0.003 


0.043 


-0.001 


± 0.001 


12.0 


± 


.4 


33.838 


20 


.473 


± 0.004 


0.044 


0.000 


± 0.001 


11.9 


± 


.4 


37.221 


20 


.682 


± 0.005 


0.040 


-0.000 


± 0.001 


14.7 


± 


.5 


40.944 


20 


.894 


± 0.006 


0.039 


-0.001 


± 0.001 


16.5 


± 


.6 


45.038 


21 


.116 


± 0.008 


0.037 


-0.002 


± 0.001 


16.6 


± 


.7 


49.542 


21 


.336 


± 0.009 


0.040 


-0.001 


± 0.001 


17.2 


±0 


.7 


54.496 


21 


.539 


±0.011 


0.041 


-0.001 


± 0.001 


18.4 


± 


.8 


59.946 


21 


.736 


±0.013 


0.044 


-0.002 


± 0.001 


19.9 


± 


.9 


65.940 


21 


.934 


±0.016 


0.040 


0.000 


± 0.001 


26.4 


± 1 


.0 


72.534 


22 


.138 


±0.019 


0.043 


-0.002 


± 0.001 


23.1 


± 1 


.1 


79.788 


22 


.328 


± 0.023 


0.048 


-0.004 


± 0.001 


23.8 


± 1 


.1 


87.766 


22 


.541 


± 0.028 


0.043 


-0.003 


± 0.001 


25.3 


± 1 


.5 


96.543 


22 


.729 


± 0.034 


0.045 


-0.005 


± 0.001 


23.9 


± 1 


,7 


106.197 


22 


.894 


± 0.039 


0.069 


-0.002 


± 0.002 


28.2 


± 1 


.5 


116.817 


23 


.031 


± 0.045 


0.081 


0.002 


± 0.002 


30.4 


± 1 


.3 


128.499 


23 


.177 


±0.051 


0.081 


0.009 


± 0.002 


37.7 


± 1 


.3 


141.349 


23 


.329 


± 0.059 


0.081 


0.000 


± 0.002 


35.4 


± 1 


.3 


155.484 


23 


490 


± 0.069 


0.097 


-0.010 


± 0.002 


28.5 


± 1 


.1 


171.032 


23 


.689 


± 0.083 


0.100 


-0.000 


± 0.002 


26.0 


± 1 


.3 


188.135 


23 


.900 


±0.102 


0.088 


-0.019 


± 0.003 


28.4 


± 1 


.7 


206.949 


24 


.076 


±0.121 


0.115 


-0.028 


± 0.003 


29.9 


± 1 


.1 


227.643 


24 


.381 


±0.164 


0.098 


-0.021 


± 0.003 


34.7 


± 1 


,7 


250.408 


24 


.611 


± 0.207 


0.113 


-0.020 


± 0.004 


38.4 


± 1 


.7 


275.449 


24 


.842 


± 0.262 


0.121 


0.002 


± 0.004 


36.4 


± 1 


.9 


302.993 


25 


.102 


± 0.346 


0.121 


-0.087 


±0.015 


45.0 


± 3 


.9 


333.293 


25 


.318 


± 0.439 


0.121 


-0.029 


± 0.007 


38.1 


± 2 


,6 


366.622 


25 


.566 


± 0.589 


0.121 


-0.001 


± 0.008 


38.1 


± 3 


.6 


403.284 


25 


.849 


± 0.854 


0.121 


0.009 


± 0.009 


38.1 


± 3 


.9 


443.613 


26 


.128 


± 1.324 


0.121 


0.027 


±0.011 


38.1 


±4 


,3 


487.974 


26 


.428 


± 2.893 


0.121 


0.069 


± 0.056 


38.1 


±4 


,8 
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Table 2. Major axis kinematics of NGC 7507. 



K 


frot 




I, 




(arcsec) 


(km s~ 1 ) 


(km s _1 ) 






-32.78 


4.4 ± 4.8 


1 7/1 1 
1Z4. 1 


± 0.0 


-0.065 + 0.027 


0.051 ± 0.017 


-28.68 


1.5 + 5.7 


1 A A 7 

144. Z 


_i_ 7 n 
± /.U 


-0.001 ± 0.028 


0.032 ± 0.020 


-25.51 


-1.5 ±5.7 


1 ^ 1 
1 JJ.l 


_i_ 7 7 

± 1 .z 


-0.001 ± 0.028 


0.025 ± 0.022 


-22.96 


-2.5 ± 6.4 


1 ^O 7 
1 Jo. / 


_i_ 7 7 
± /.Z 


0.016 ± 0.028 


-0.003 ± 0.025 


-20.78 


-11.8 ±5.7 


1 A/1 A 
104.4 


± o.y 


0.043 ± 0.026 


0.011 ±0.026 


-18.88 


3.5 ± 5.8 


1 A/1 1 
104.J 


_i_ O A 
± 0.0 


-0.007 ± 0.027 


0.057 ± 0.025 


-17.28 


-9.4 ± 6.3 


1 QCi A 


_i_ O A 
± 0.0 


0.007 ± 0.024 


0.047 ± 0.025 


-15.90 


4 ±5.8 


1 HZ. A 
1 /J.O 


_i_ 7 n 
± 1 .y 


-0.033 ± 0.025 


0.048 ± 0.022 


-14.66 


-5.8 ± 6.7 


1 QQ 


-i- 7 Q 


-0.002 ± 0.024 


0.018 ± 0.025 


-13.57 


9 ±6.9 


1 07 O 

10 1 .y 


-i- 7 Q 


-0.010 ± 0.024 


0.039 ± 0.024 


-12.62 


-3.3 ± 6.0 


1 Q1 7 


-i- A 
± 0.0 


0.012 ± 0.023 


0.046 ± 0.023 


-11.82 


1.2 ± 6.5 


1 OA fl 


_i_ Q 7 
± O.Z 


0.012 ± 0.025 


0.015 ± 0.023 


-11.09 


-4.4 ± 7.6 


1 Q7 1 

iyZ. 1 


-u 1 
± 0. 1 


-0.007 ± 0.024 


0.021 ± 0.026 


-10.43 


-3.1 ± 6.9 


1 QA "X 


± O.J 


-0.015 ± 0.024 


0.025 ± 0.024 


-9.85 


-3.4 ±6.1 


1 07 1 
10 / . 1 


-i-Q7 
± O.Z 


-0.007 ± 0.023 


0.036 ± 0.023 


-9.26 


-8.8 ±5.8 


1 m 7 
iy 1.1 


_i_ 7 n 
± i.y 


-0.000 ± 0.022 


0.018 ± 0.023 


-8.75 


-2.1 ± 7.0 


1 7Q 7 
1 ly.Z 


± O. / 


-0.015 ± 0.026 


0.053 ± 0.023 


-8.31 


1.8 ±6.7 


1 Q 

iyj.y 


_i_ 7 7 
±1.1 


0.008 ± 0.024 


0.015 ± 0.022 


-7.88 


-0.1 ± 5.9 


7 1 C\ 1 

Z1U. 1 


_i_ 7 A 
± /.0 


-0.004 ± 0.022 


0.041 ± 0.024 


-7.44 


6.3 ± 6.4 


1 CK A 

IVj.o 


± O.U 


0.014 ± 0.021 


0.066 ± 0.024 


-7.00 


-7.2 ± 5.4 


Zl 1.0 


_i_ O 1 
± 0. 1 


-0.038 ± 0.020 


0.034 ± 0.025 


-6.64 


0.2 ± 7.0 


7ha a 
ZU0.4 


_i_ Q 7 

± o.Z 


-0.005 ± 0.023 


0.043 ± 0.026 


-6.05 


-6.9 ±6.1 


1 fl/1 A 

1V4.4 


_i_ 7 Q 
± /.0 


-0.006 ± 0.022 


0.088 ± 0.022 


-5.47 


1.4 ±6.4 


7flQ A 

Zvy.O 


-u 7 1 
± /.I 


0.043 ± 0.019 


0.034 ± 0.022 


-5.18 


5.7 ±6.1 


Z1Z.4 


_i_ 7 7 
±1.1 


-0.025 ±0.019 


0.076 ± 0.020 


-4.59 


7.8 ± 5.5 


7 i n 7 

Z1U.Z 


-L- 7 1 
± /.I 


0.006 + 0.018 


0.053 ± 0.019 


-4.01 


-3 ±5.1 


71 A Q 
214.5* 


_i_ A 7 

± 0. / 


0.038 + 0.015 


0.064 ± 0.019 


-3.50 


-0.1 ± 6.5 


Z1Z.0 


± o.U 


-0.012 ± 0.017 


0.065 ± 0.022 


-3.06 


-0.2 ± 5.9 


770 Q 

ZZo.y 


_i_ 7 ^ 
± /.J 


0.052 ± 0.016 


0.041 ± 0.021 


-2.04 


-1.7 ±4.8 


ZlO.O 


+ J.Z 


0.005 ± 0.014 


0.062 ± 0.015 


0.00 


-0.2 ± 3.5 


7/1 C A 

Z40.4 


± 4.U 


0.026 ± 0.009 


0.055 ± 0.011 


1.02 


2± 3.8 


Zjj.Z 


± 4.4 


-0.010 ± 0.012 


0.051 ± 0.013 


2.04 


7.8 ±5.3 


ZZv. 1 


_i_ A 7 
± O.Z 


-0.009 ± 0.017 


0.037 ± 0.015 


3.07 


-0.2 ± 6.2 


ZZvA 


± O.J 


0.012 ± 0.022 


0.080 ± 0.022 


3.50 


7.1 ± 7.5 


1 C 1 

Zlj. 1 


± 0.4 


0.011 + 0.025 


0.065 ± 0.021 


4.16 


3.2 ± 6.5 


1 f\ 1 

Z1U. / 


_i_ A O 

± o.y 


-0.001 ± 0.022 


0.056 ± 0.020 


4.45 


0.9 ± 6.0 


1 1 1 A 


_i_7 
± /.0 


-0.006 ± 0.023 


0.036 ± 0.021 


5.03 


-1.5 ±7.0 


ZUU. / 


± O.J 


0.000 ± 0.025 


0.047 ± 0.026 


5.62 


1.0 ±6.8 


ZUj.U 


_i_0 7 

± o.Z 


0.004 ± 0.025 


0.045 ± 0.023 


5.91 


4.4 ± 7.4 


OHQ 1 

ZUo\ 1 


± O.J 


-0.041 ± 0.026 


0.060 ± 0.023 


6.50 


-5.5 ±7.8 


ZUo.i5 


_i_ O 7 
+ 0. / 


-0.011 ±0.028 


0.007 ± 0.026 


6.86 


7.9 ± 6.7 


1 no 1 
iVo. 1 


_i_ 7 O 
± /.0 


-0.005 ± 0.024 


0.057 ± 0.022 


7.73 


-4.9 ± 6.7 


1 nQ 1 

iyj. i 


_i_ Q 7 
± O.Z 


0.007 ± 0.027 


0.023 ± 0.023 


8.17 


-12.8 ±7.2 


1 C7 ^1 
lo / . J 


± y.u 


-0.008 ± 0.028 


0.055 ± 0.025 


8.61 


3.2 ±7.3 


1 m a 
iy /.o 


_i_ A 
± 0.0 


-0.020 ± 0.028 


0.049 ± 0.022 


9.12 


-7.5 ± 6.6 


loo. J 


_i_ O 1 
+ 0.1 


-0.028 ± 0.027 


0.047 ± 0.023 


9.70 


0.5 ± 7.2 


1 OQ 7 


± O.J 


-0.005 ± 0.026 


0.033 ± 0.023 


10.29 


-13.4 ±7.9 


1 ^ 


_i_ o n 


0.032 ± 0.027 


0.032 ± 0.023 


10.94 


-6.4 ± 7.0 


loi.y 


± O.J 


0.006 ± 0.028 


0.059 ± 0.022 


11.67 


1.5 ± 8.2 


1 oO.j 


± 0.0 


0.021 ± 0.029 


0.056 ± 0.024 


12.47 


2.6 ± 7.2 


1 1 

iyj.i 


_i_ o ^ 
± O.J 


0.005 ± 0.028 


0.058 ± 0.021 


13.35 


-5.0 ± 7.6 


1 7/1 c 
1 /4.0 


-i- Q 7 
± y.Z 


-0.033 ± 0.029 


0.071 ± 0.028 


14.30 


1.9 ±7.9 


1 OJ.J 


± y.J 


0.013 ± 0.029 


0.036 ± 0.026 


15.39 


1.4 ±7.8 


1 Q 1 7 
IV 1 .Z 


-L- 

± 0.0 


0.004 ± 0.029 


0.010 ± 0.027 


16.63 


-5.9 ±7.8 


1 on a 
loU.O 


± 0.4 


0.037 ± 0.030 


0.046 ± 0.023 


18.01 


-8.0 ± 7.2 


188.6 


±8.5 


-0.006 ± 0.030 


0.001 ± 0.026 


19.62 


1.3 ± 8.0 


177.5 


±7.9 


-0.042 ±0.031 


-0.016 ± 0.025 


21.44 


-13 ±7.4 


172.8 


+ 8.1 


-0.003 ±0.031 


0.029 ± 0.024 


23.48 


-15.4 ±7.3 


155.6 


+ 8.3 


0.001 ± 0.031 


0.022 ± 0.024 


25.88 


-6.2 ±7.8 


148.4 


±7.8 


-0.030 ± 0.032 


0.053 ± 0.022 


28.78 


-2.8 ± 6.8 


139.2 


±7.8 


-0.034 ± 0.029 


0.083 ± 0.019 


32.40 


-21.1 ± 6.9 


133.4 


+ 8.1 


0.068 ± 0.029 


0.077 ± 0.020 


37.30 


-11.4 + 7.0 


146.9 


± 7.4 


0.024 ± 0.031 


-0.004 ± 0.025 


44.70 


-12 ± 6.8 


143.4 


±7.3 


-0.001 ± 0.027 


0.030 ± 0.017 


58.08 


-19.3 ±6.6 


116.7 


±8.5 


0.079 ± 0.017 


-0.012 ± 0.012 


97.33 


-68.0 ± 8.1 


109.6+18.3 


-0.009 ± 0.017 


-0.025 ± 0.018 
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Table 3. Minor axis kinematics of NGC 7507. 





^rot 


0~y 






tl4 


yen ^ scl f 


(km s -1 ) 


(km s~') 








-81.02 


-21.64 ± 26.2 


133.5 ± 32.8 


-0.021 ± 0.007 


-0.004 ± 0.006 


-46.29 


25.06 ± 10.2 


156.8 ± 11.5 


-0.077 ± 0.019 


U.UUj 


, Af>1Q 

± u.uiy 


-34.90 


22.56 ± 8.5 


163.3 ± 10.1 


0.008 ± 0.023 


n no/i 
U.UV4 


± U.Ulo 


-27.88 


42.26 ± 8.0 


159.9 ± 9.7 


0.018 ± 0.023 


n 1 oo 


+ U.U1 J 


-16.66 


34.56 ± 7.2 


184.5 ± 9.5 


0.017 ± 0.026 


u.uoy 


+ U.UZ4 


-14.34 


28.96 ± 7.2 


173.6 ± 7.5 


-0.016 ± 0.026 


n 

U.Ujd 


± u.uzz 


-12.53 


32.26 ± 7.1 


193.2 ± 9.4 


-0.004 ± 0.026 


n noQ 
U.Uoj 


± U.UZz 


-11.07 


33.36 ± 7.5 


196.5 ± 9.6 


-0.064 ± 0.026 


n an 
U.Ulz 


± U.UZo 


-9.84 


29.56 ± 7.2 


185.9 ± 8.6 


-0.046 ± 0.026 


U.U4/ 


± U.UZj 


-8.74 


25.46 ± 6.3 


184.5 ± 9.1 


0.040 ± 0.025 


n 

U.Ujj 


_i_ n fin 
± U.UZ / 


-7.80 


19.86 ± 5.8 


177.5 ± 8.4 


0.021 ± 0.025 


fl f\HQ 

U.U/o 


± U.UZ4 


-7.00 


33.46 ± 6.1 


188.8 ± 7.9 


-0.036 ± 0.024 


n n^A 
U.Ujo 


i n noo 


-6.35 


27.86 ± 7.2 


194.4 ± 8.3 


-0.040 ± 0.026 


n no< 
U.Uzj 


, f\ no^ 
± U.Uzj 


-5.76 


21.96 ± 6.3 


201.5 ± 8.4 


0.007 ± 0.025 


n ni 1 
U.Uj 1 


_i_ n no i 


-5.18 


18.76 ± 5.8 


198.5 ± 7.5 


-0.009 ± 0.024 


U.U4o 


_i_ n noo 
± U.UZz 


-4.67 


15.76 ± 6.4 


205.7 ± 8.9 


0.011 ± 0.024 


n n^c 
U.Ujo 


± U.Uzo 


-4.23 


16.76 ± 5.9 


211.4 ± 8.1 


0.020 ± 0.022 


n r\ir\ 


_i_ n no i 
± U.Uzj 


-3.79 


20.16 ± 5.0 


201.2 ± 8.0 


0.029 ± 0.022 


U.UoU 


i n noo 
± U.UzZ 


-3.43 


14.56 ± 6.0 


204.1 ± 9.2 


-0.007 ± 0.024 


n ci&.i 
U.Uoo 


_i_ n no^ 
± U.Uzj 


-3.14 


16.56 ± 6.1 


210.7 ± 8.4 


0.005 ± 0.022 


U.U40 


i n noQ 
± U.Uzj 


-2.04 


1.26 ± 6.2 


217.9 ± 8.8 


0.035 ± 0.023 


U.Ujj 


, r\ noo 
± U.UZZ 


-1.02 


7.76 ±4.8 


227.4 ± 6.7 


0.022 ± 0.017 


U.Uo4 


_i_ n m 6. 
± U.Ulo 


0.00 


8.36 ± 4.0 


245.4 ± 5.3 


0.029 ± 0.012 


C\ (Y7Q 

U.U/o 


_i_ n n 1 1 
± U.U1 5 


1.02 


-11.04 ± 4.9 


224.5 ± 5.3 


0.034 ± 0.016 




_i_ n m n 


2.11 


-10.44 ± 4.3 


212.1 ± 5.9 


0.019 ± 0.015 


U.Uu/ 


_i_ n non 
± U.UZU 


3.28 


-20.44 ± 6.0 


207.2 ± 7.1 


-0.003 ± 0.020 


u.u iy 


_i_ n no a 
+ U.UZ4 


3.64 


-26.94 ± 5.1 


209.8 ± 6.9 


0.029 ± 0.017 


U.U42 


_i_ n non 
± U.UZU 


4.08 


-23.24 ± 6.4 


212.7 ± 6.6 


0.035 ± 0.019 


U.UJU 


i n no i 
± U.Uzl 


4.52 


-23.14 ± 5.9 


203.1 ± 6.6 


-0.001 ± 0.020 


r\ m 1 


, f\ noo 
± U.UZZ 


5.03 


-33.94 ± 5.2 


202.8 ± 6.4 


0.013 ± 0.019 


U.UIJ 


i n noQ 
± U.Uzj 


5.62 


-22.94 ± 6.9 


189.9 ± 7.7 


-0.043 ± 0.023 


U.U3V 


, f\ no^; 
± U.Uzo 


6.27 


-32.24 ± 6.5 


183.9 ± 7.5 


0.023 ± 0.023 


U.Uoz 


, f\ noQ 
± U.Uzj 


7.00 


-30.14 ± 6.2 


204.0 ± 7.4 


-0.020 ± 0.022 


0.021 


± 0.027 


7.80 


-36.94 ± 6.6 


173.3 ± 7.3 


0.014 ± 0.023 


0.104 


±0.021 


8.74 


-35.04 ± 6.2 


187.7 ± 6.4 


0.069 ± 0.023 


0.023 


± 0.025 


9.84 


-54.44 ± 7.4 


186.8 ± 7.0 


0.079 ± 0.023 


0.064 


± 0.024 


11 i/i 

11.14 


-51.14 ± 1.0 


1 /9.9 ± I.Z 


U.04/ ± I). 1)24 


0.069 


± 0.024 


12.81 


-43.04 ± 6.6 


184.9 ± 7.9 


0.018 ± 0.024 


0.089 


± 0.025 


14.92 


-54.44 ± 6.6 


174.4 ± 7.3 


0.035 ± 0.025 


0.056 


±0.021 


17.53 


-61.14 ±6.5 


173.2 ±7.7 


0.033 ± 0.025 


0.051 


± 0.025 


20.87 


-55.14 ±7.5 


179.3 ± 6.9 


0.061 ± 0.025 


0.015 


± 0.024 


32.23 


-54.74 ± 8.3 


158.2 ± 8.3 


0.047 ± 0.024 


0.031 


± 0.023 


42.16 


-48.24 ± 7.5 


122.1 ± 11.9 


0.027 ± 0.014 


0.034 


±0.013 


61.48 


-8.44 ±21.1 


135.9 ± 27.4 


-0.021 ± 0.009 


-0.013 


± 0.007 



